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Abstract. We discuss the propagation of hydrogen atoms in static electric and magnetic fields in a lon- 
gitudinal atomic beam spin echo (lABSE) apparatus. There the atoms acquire geometric (Berry) phases 
that exhibit a new manifestation of parity-(P-)violation in atomic physics. We provide analytical as well 
as numerical calculations of the behaviour of the metastable 2S states of hydrogen. The conditions for 
electromagnetic field configurations that allow for adiabatic evolution of the relevant atomic states are 
investigated. Our results provide the theoretical basis for the discussion of possible measurements of P- 
violating geometric phases in lABSE experiments. 
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1 introduction 

In this paper we continue our theoretical studies towards 
the realisation of an experiment to measure P-violating 
effects using the methods available with an atomic beam 
interferometer as described in [1]. Starting with the pa- 
per [2] many theoretical and experimental studies of P- 
violation due to neutral current exchange were made for 
heavy atoms. For reviews see [3,4]. The most precise re- 
sults so far are from experiments with Cs [5,6]. The exper- 
imental observation of P-violation in the lightest atoms, 
hydrogen (H) and deuterium (D), still is an open prob- 
lem [7] . For the early experimental efforts in this direction 
see [8-10]. In [11, 12] we have studied parity violating geo- 
metric phases. These are interesting quantum mechanical 
phenomena by themselves and we think it is worthwhile 
to try measuring them. In the present paper we shall de- 
velop the necessary theoretical tools for a description of 
the behaviour of an atom in the interferometer described 
in [1]. The aim is to calculate the relevant phases which 
atoms pick up when travelling through the interferome- 
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ter and to construct observables sensitive to P-violating 
geometric phases. We shall focus on the hydrogen states 
with principal quantum number n = 2, in particular on 
the metastable 2S states. Our paper is organised as fol- 
lows. In Section 2 we set up the basic equations governing 
the motion of the atom in the n = 2 states as it trav- 
els through the interferometer. In Section 3 we study the 
equation for a single substate. Section 4 deals with two 
metastable 2S states and with their interference effects. 
In Section 5 we draw our conclusions. A concrete pro- 
posal for a measurement of P-violating geometric phases 
will be given in a follow up paper. Our notation is ex- 
plained in Appendix A which collects also many useful 
formulae. Furthermore, we discuss in Appendix A the con- 
ditions for adiabatic evolution of the atomic states for the 
case considered. Appendix B contains the details of our 
calculations concerning the motion of a single 2S state in 
the interferometer. In Appendix C we discuss in detail the 
interference effects of two states. Throughout our work we 
use natural units h = c = I if not indicated otherwise. 



2 A hydrogen atom in the interferometer 

We consider in this section the case of a hydrogen atom 
being in a superposition of n = 2 internal states and trav- 
elling through an interferometer of the type presented in 
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[1]. The latter is, for a theorist, a device producing static 
electric and magnetic fields acting on the atom. Schemati- 
cally this is shown in Figure 1. The magnetic field strength 
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Fig. 1. Schematic view of the longitudinal spin echo interfer- 
ometer. The coordinate system used is indicated. The coor- 
dinate along the beam direction is denoted by z. The atom 
is prepared in a wave packet around zo and analysed around 
Za- In the example discussed in Section 4 we start with a su- 
perposition \iIj{zo)) of two states, labeled with |9) and |11). 
After passing the magnetic and electric fields B{z) and £{z) 
the wavefunction is projected onto an analysing state \'4>{za)), 
for example again onto a superposition of the states |9) and 
111). 



3 and the electric field strength £ are supposed to depend 
only on the z coordinate and to be nonzero only between 
Zo and Za- Here zq marks the place at the beginning of 
the interferometer around which, in a field-free region, the 
wave packet of the atom is prepared. Similarly, Za at the 
end of the interferometer marks the place where, again in 
a field-free region, the internal state of the atom is anal- 
ysed. 

The hydrogen atom in a superposition of states with 
principal quantum number n = 2 is described by a 16- 
component spinor 



16 



|^(x,t)) = ^/„(x,t)|a;0) . 



(1) 



Here |a; 0), a = 1, . . . , 16, arc the basis states for 3 
and 5 = 0, that is the 2Si/2, 2Pi/2 and 2P 



3/2 



states 



with the numbering as given in Table 2 of Appendix A. 
The Schrodingcr equation describing the undccaycd n = 2 
states reads 



1 

2m 



A + .ja,z) + 



k'(x,t)). (2) 



Here m is the atom's mass and .jMi^z) is the non-hermitian 
energy or mass matrix as described in Section 2.2 of [12]. 
There the center of the real parts of the energy levels of 
the free 2Pi/2 states was chosen as the (arbitrary) zero 
for the energy scale. In the following it is convenient to 
choose this zero point of the energy differently. Therefore, 
we introduce in (2) a real constant which we can choose 
freely. The mass matrix ^^(z) is given by 



(3) 

(4) 

Here we use the same notation as in (22)-(24) of [12]. 
The mass matrix for zero external fields and without P- 
violating contributions is denoted by The P-violating 



parts of the mass matrix are J^^py , z = 1,2. The two 
P-violation parameters bi quantifying the nuclear spin de- 
pendent and independent parts, respectively, are defined 
in Appendix A. The matrices D_ and are the electric and 
magnetic dipole-moment operators for the n = 2 states. 
The explicit form of all these matrices is given in Appendix 
D. 

Let the momentum of the atom be 



<7t = 



1y 



(5) 



The problem is to discuss the solution of (2) for the case 
that the atom travels essentially in z direction where for 
all z the inequality 



2m 



holds. Here ||i^(z)|| is a suitable norm of ^{z) 

1/2 



||^(z)|| = tr[^\z)M.{z) 



(6) 



(7) 



Similarly, the kinetic energy of transverse motion is always 
supposed to be much smaller than that of longitudinal 
motion 



2m 2m 



(8) 



Furthermore we suppose ^M_(z) to vary only smoothly with 
z, such that the conditions for adiabatic evolution of the 
relevant states are fulfilled. This is specified for a concrete 
example below, and we discuss the adiabaticity criterion 
in detail in Appendix A. For z around zq and around Zq 
(see Figure 1) we suppose absence of external fields, thus 



£{z) = 0, S(z) = 



(9) 



for 



\z~za\<8l and \z - Za\ < 51 ■ (10) 
Here 81 is supposed to be very small compared to Za — 20, 
< (5/ < Za - Zo . (11) 

We have then 

^(z) = ^(zo) = MiZa) = ^ol-^i, ^2) (12) 

for \z — zq\ < bl and |z — Za| < bl. 

The next step is to introduce the right and left eigen- 
vectors of the non-hermitian matrix ^M_(z). We assume 
here that ;M_{z) can be diagonalised for all z. For each 
z we have then 16 linearly independent right and. analo- 
gously, left eigenvectors 



^(z)|a(z)) = S,(z)|a(z)) , 
[^z)\M.z)^{fAz)\E^{z) , 



,16 . 



(13) 
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The eigenvalues Ea{z) are in general complex. We nor- 
malise the eigenvectors such that 



{a{z)\[3{z)) ^ Sa(3 



and 



{a{z)\a{z)) = 1 , 

(no summation over a) . 

We can now expand \'ip{x,t)) in terms of 

16 

|V'(x,t)) =^7/.a(x,t)|a(z)) . 
From (2) we get then 

16 



(14) 
(15) 

(16) 



E 

13=1 



-^A^p{x,t) + {Ep{z) + So) V/3(a;,t) 
Zm 



1 (di^a 



l/3(^)) 



3z 



-^#(a.,t)^|/3(.)) 







(17) 



Multiplying with (a(z)| from the left gives 16 coupled 
equations 

^ ZiVa(a;, t) + (-B„(z) + So) ^a(a;, t) - i^(a^, 



2m 

02 m az 

13=1 ^ 



d 



1 



(a{z)\^\(3iz))—^pix,t)\^0 



2m 



(a = l,...,16) 



(18) 



This gives a kinetic energy of longitudinal motion 
= 6.5 X 10^2 eV 

2m 



(21) 



and a time of flight of the atom through the interferometer 
of 



1.43 ms . 



(22) 



The lifetimes of the free 2S and 2P states are [13,14] 
Ts = Fg^ = 0.1216 s , 



TP = Fp' 



1.596 X 10"^ s 



(23) 



In an electric field the 2S states obtain a 2P admixture 
and their lifetimes decrease depending on the strength of 
S, see Figure 2. To have an appreciable flux of metastable 



0.1 



0.01 




In the following we shall consider (18) in the adiabatic 
limit. That is, we want to choose conditions where the 
evolution of the metastable 2S states (a = 9, . . . , 12) de- 
couples from that of the 2P states, see [11,12]. This gives 
restrictions on the allowed rate of change of£{z) and B{z) 
with z. In order to get physical insight about these restric- 
tions we shall now analyse them for a concrete case. Of 
course, the methods and formulas given below are easily 
adapted to other cases. 

In the following example, the interferometer has a length 

of 



0.001 



0.1 0.2 0.3 0.4 0.5 0.6 

\£\ [V/cm] 

Fig. 2. Lifetimes of the mixed 2Si/2 states of hydrogen in the 
external electric field £ = £es as function of \£\. The number- 
ing of the states (a = 9, . . . , 12) is explained in Table 2, see 
Appendix A. 



2S states at the end of the interferometer one should in 
view of (22) require, at least on average. 



Ts{S) > 1 ms 



(24) 



From Figure 2 we see that this limits the allowed electric 
fields on average to 



I 



zo « 5 m 



(19) 



\£{z)\ay < 0.60 V/cm 



(25) 



The typical longitudinal velocity w^, wave number Qz and 
de Broglie wavelength A of the H-atom in the interferom- 
eter at room temperature, 300 K, are 



Of course, locally the fields can be stronger, but in practice 
a reasonable limit seems to be 



3500 m/s , 



|f(^)|_,<10V/cm 



(26) 



5.6 X 10^° m^^ 



1.1 X 10"^"m . 



(20) 



Estimates for the allowed rate of change of S{z) and 
13{z) compatible with adiabatic evolution of the atomic 
states are discussed in Appendix A based on the results 
of [12]. We find that the evolution of the 2S states in the 
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interferometer will decouple from that of the 2P states if 
the typical length scale Az for a significant variation of 
the electric field S{z) is Az > 50 /im. For the magnetic 
field B{z) the length scale of significant variation should 
satisfy Az > 10 /um, see (A.15)-(A.22) of Appendix A. 

The occurrence of transitions amongst the individual 
2S hyperfine states a = 9, . . . , 12 in the interferometer 
depends very much on the precise setup of magnetic and 
electric fields since for B{z) = and £{z) = the states 
a = 9, 10 and 11, that is, the states 2Si/2 with F = 1 and 

= ±1,0, are degenerate. The Brcit-Rabi diagram for 
the states a = 9, . . . , 12 is shown in Figure 3. For concrete 
cases one must check the adiabaticity conditions explicitly. 
Typically the requirement is that the energy differences 
AE in the Breit-Rabi diagram for the B-fields considered 
correspond to frequencies larger than the frequencies 1/At 
of variation of the £- and B- fields. Setting Az = VzAt for 
the corresponding variation in z we get 



AE > 



Az 



Az > 



AE/h 



(27) 



In this way we get, for example, for B = Be^ with B = 
0.5 mT for the = 1 states from Figure 3 AE 0.05 ficV, 
AE/h w 12 MHz, Az > 300 ^m. Of course, for other re- 
gions of the Breit-Rabi diagram one can make correspond- 
ing estimates. 
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Fig. 3. Breit-Rabi diagram for the 2S states of hydrogen. The 
magnetic field is B = Bes. 



Thus, we see that for the above example the conditions 
for adiabatic evolution of the 2S states lead to experimen- 
tally realisable restrictions for the electric and magnetic 
fields. With this we end our discussion of the example 
here, but we shall return to it in later sections. 

Having discussed exemplarily the constraints on the 
electromagnetic field configuration for the adiabatic ap- 
proximation to hold, we continue with the general inves- 
tigation of (18). In the following we suppose that the con- 
ditions for adiabatic evolution of the metastable 2S states 
are fulfilled. Then, the evolution of these 2S states is gov- 
erned by four decoupled equations obtained from (18) for 
a = 9, . . . , 12 and keeping there in the sum over /3 only 



the terms with (3 = a. In this way we get as the basic 
system of equations for the following discussions 

- J-A^p^{x,t) - -{c^)\^\aiz))^ix,t) 
Zm m oz oz 



Eciz) + Eo - :^{a{z)\-^\a{z)) 



2m 



d 

i^^a(a;,t) = , 



(a = 9,...,12) . 



(28) 



To write (28) in a convenient way we introduce the diag- 
onal geometric phases (see (45) of [12]) 



d 



-faa{z) = i / dz' {a{z')\ — \a{z')) , (a = 9 



12) 



(29) 



Note that for metastable states jaaiz) is, in general, com- 
plex. We split off from the wave function a factor corre- 
sponding to this geometric phase and define 



Tpaix.t) = exp[-i7aa(2)]V'a(a;,t) 

For Tp we find from (28) the equation 



(30) 



AipJx^t) + 2mi^$^(x,t) - 2mVJz)ipJx,t) = (31) 
ot 



where we define 



V^iz) ^ E^iz) + Eo - — 
Zm 



+ 1- 



dz 



dz^ 



— 52 
+ Hz)\-^\aiz)) 

E^{z) + Eq 

'd-fc.c.{z)Y d{^)\d\a{z)) 



1 
2m 



dz 



dz dz 



(32) 



In the following we set 



vM)^v^iz)~-r^{z) 

y«(z) = ReV„(z) , 
r,(z)--2ImVa(z) . 



(33) 



The effective potential Vq,(z) is invariant under local phase 
transformations of the eigenstates |q!(2^)) according to 

|a(z))-^|a(z))' = e"'=(^)|a(z)) , 
(a(^)|^(a(I)|' = e-''°(^)(a(I)| (34) 

where rja{z) has to be real in order to respect (15), see 
Appendix B. 

The solutions of (31) will be discussed in the next sec- 
tions using methods from wavelet theory, see for instance 
[15]. 
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3 Solution for one component 

We now study the solution of the effective equation (31) 
for the amphtudes of the metastable states for one com- 
ponent a. Since we keep a fixed in this section we omit 
this index from our quantities in the following. Thus, we 
drop the index a in (30) 

{j:{x, t) = exp[-i7(z)]V'(a;, t) , (35) 

and the equation to be studied for ip reads 



^2 Q 

^ + 2im--2my(z) 



+ imr(z) + Z\TU(a;,t) = (36) 



effect in the wave function when k{z) changes as well as the 
decay of the metastable atom, we rather choose a complex 



[z,t) 



where 



2m 



dz' fc(z') 



iln^(^) 
2 



1 

+ 2 



dz 



k{u) 
,mr{z') 
k(z') 



u{z,t)=Z{C-T) 



Z{^{z)--t 



u—u{z,t) 

(43) 



(44) 



with 



XT 

z 



At = 



(37) 



Let fc, of order g^, be a typical scale for the momentum 
component in z direction. We change variables and set 



k 

—t 

m 



dz'- 

Jk'^-2mV{z') 



(38) 



From (38) we get C as function of z. For this and for the 
inverse function we write 



C = 3(^) : 

z = Z{C) 



(39) 



Furthermore, we split off from the wave function ^{x, t) 
a convenient factor exp[i(/)(z, t)] and denote by A the re- 
maining amplitude which is the envelope function for the 
wave packet: 



i,{x,t) = e'^^'^'^A{xTX.r) 



(40) 



In the following we could choose to be the WKB phase 
factor 



4>WKS,{z, t) 



t+ / dz'fc(z') 



Here, the local momentum k{z) is defined as 



k{z) = Jk'^ - 2mV{z) 



(41) 



(42) 



Then exp[i(j>^/vKB{z,t)] is a pure phase factor since (/)wkb 
is real. However, taking into account a simple kinematic 



It is now a simple exercise to transform (36) to the 
new variables xttC^t and to derive the equation for the 
envelope function A{xt,Cit) of (40). For the details see 
Appendix B. The result can be written in the form 



d_ 
dr 



d_ 



A{xt,C,t) = {LA){xt,C,t) 



(45) 



where L is an operator given in detail in (B.21)-(B.24). 
The initial condition for A is 



(46) 



with some given function ip determined by the initial ex- 
perimental conditions. With this the solution of (45) is 
easily obtained by transforming (45) into an integral equa- 
tion. Here we use a Green's function method, see Ap- 
pendix B. The final result is an expansion of A, see (B.39), 



^(a:T,C,T) = ^^("na;T,CT) (47) 

n=0 

where successive terms are suppressed by higher and higher 
powers of 1/fc 



^(")(a3T,C,r)=0(lA") 



(48) 



For large k the zero order term should be a good approx- 
imation. It reads, see (B.40), 



A^"\xT,C,T)^ip{xT,C-r) 



(49) 



The higher order terms in (47) are discussed in Appendix 
B. The zero order term (49) describes the motion of the 
wave packet without longitudinal and transverse disper- 
sion. The envelope function is just progressing with con- 
stant speed in the ( coordinate, where C ~ 5{z), see (38) 
and (39). 

To discuss the physics of the solution (49) let us sup- 
pose that for t = the envelope function (p{xT,3{z)) is 
sharply peaked at xt = and z = zq corresponding to 
C = 0. Then, according to (49), the envelope function for 
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times t > stays equally peaked at xt = but the peak the roadmap plot is shown for the states a = 9 and a = 11 



in z moves according to z(t) where 



3(zW)--i = 3(^o) = o. 

m 

Inserting here 3(0 from (38), (39) we get 



=(t) 



dz' 



fc2 - 2mV{z') 



(50) 



(51) 



which gives 



-mz'^{t) + Viz{t)) 



2m 



(52) 



The motion of the peak in z corresponds to that of a 
classical particle in the real part of the potential where 
the total energy is W /2m. 

Now we can make contact with the roadmap descrip- 
tion ( "Fahrplanmodell" ) , introduced in [16] and in a manu- 
script currently prepared for publication [17]. We consider 
wave packets starting at t = in the field free region at zq, 
see Figure 1. We define the arrival time t{z) and the corre- 
sponding reduced "time" t{z) by reversing the functional 
dependence in (51) 



t{z) 



fc2 - 2mV{z') 



r{z) = -t{z) ^ 5{z) . 
m 



(53) 



Note that the reduced "time" r has the dimension of 
length. We compare these arrival times with those of a 
free wave packet where V{z') = 0. From (32) and (33) we 
see that this corresponds to a free atom where "f{z) = 
and the energy level E{z') satisfies 



ReE{z) + Eo = 



(54) 



That is. the choice of Eq selects the free state to which one 
compares, where the arrival time t and reduced "time" r 
are 



m 
X 
k 



{z - ^o) , 

tfroc(z) = Z - Zo . 



(55) 



The difference of reduced arrival times 

At{z) = t{z) - rfi.co(z) 
= 3(^) -{z- Zq) 



(56) 



describes the delay [At > 0) or advance {At < 0) of 
the arrival of the wave packet peak at coordinate z in the 
interferometer with respect to that of the free atom. 

As a concrete example we consider again the interfer- 
ometer with the parameters specified in (19) ff. In Figure 5 



in a purely longitudinal magnetic field 
B{z)=B{z)e3 , 

Biz) = 10 (e-6(--i-25)' - e-6(-3.75)= j , (57) 
depicted in Figure 4, where z is taken in meters. We com- 

12 



42 




z [m\ 

Fig. 4. Example of a magnetic field configuration. The trans- 
verse components are chosen as Bi = B2 = and the longitu- 
dinal component B{z) is given in (57). 

pare the evolution of these states with that of the free 2S, 
F = 1, states (a — 9, 10 and 11). That is, we set 

^0 = - ^ , (58) 

see Table 3 in Appendix D. Here A denotes the ground 
state hyperfine splitting energy. The potentials for the 
states a = 9 and a = 11 are then proportional to the 
third component of the magnetic field in (57): 

5MB , 



Vgiz) 

Vniz) 



2 



Biz) , 
Biz) 



(59) 



with the Bohr magneton fiB and the electron's Lande fac- 
tor g, see Table 5 in Appendix D. We see from Figure 5 
that in the magnetic field (57) the state a = 9 first gets 
a delay, a = 11 an advance compared to the free state. 
For z around 3.75 m where the magnetic field is reversed 
the state a = 9 makes up for the delay, the state a = 11 
for its advance. Finally at z = 5 m they meet again. This 
gives a nice visualisation of the motion of the wave packet 
in the interferometer. Clearly, the complete solution con- 
tains much more information, as discussed below. 

Returning now to the general discussion and putting 
everything together we find from (35), (40), (43) and (49) 
the lowest order solution for the wave function as 



ipix, t; k) = exp[i7(z)]'(/;(a;, t) 

= exp[i0(z, t) + i7(z)]^(") ixT, C, t) 
= exp[i0(z, t) + i7(z)].^(a;T, ( ~ t) . 



(60) 
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Fig. 5. Example of a roadmap plot for the states a = 9 and 
Q = 11. The corresponding potentials are given in (59). 



For purposes which will become clear in Section 4 we in- 
dicate here explicitly also the dependence of ^ on k. To 
write this in a physically transparent form we note that 
(see (42)) 



k{z) 



2mV{z) 
'k{z) + k ' 



dz'k{z') ^k{z- zo)+ dz'(k{z')~k) 



U \ r A /2mF(£0 



(61) 



We define now 

[/(z, m; fc) = exp 



-ilnMf)_i 
2 kiu] 2 



dz 



,2mV{z') 
k{z') + k 
^,mr{z') 



(62) 



It is easy to see that U corresponds to the evolution of the 
wave function of an atom at rest with an effective complex 
mass 



i dV{z) 
2k{z) dz 



2 m dz 



(63) 



We get 



k[z) dU{z,u;k) 
m dz 
U{z,u;k)\ 



-i^„g(z)U(z, u; k) 



(64) 



The solution (60) can now be written as 



-i— t + ik{z - zo) 
Zm 



ip{x, t; k) = exp 

X U{z,u;k)\^^^^^^_^^x{xr,C, ~ t) 



(65) 



where we define 

x{xtX-t) = 

exp 



dz'— — + i7(u) 



k{z') + k 
X '^{xtX-t) . 



(66) 



For T = we have u = Z{C,) = z and U{z, u;k) = \. Thus, 

V'(a;,t;fc)|j^o = exp[ifc(z - Zo)] x{xt,0\c=^{z) (67) 

contains a plane wave factor and x(a;T,C)i the envelope 
amplitude at i = 0. The interpretation of the result (65) 
is then as follows. The wave function ip{x,t;k) contains 
a plane wave factor and a factor U{z,u) corresponding 
to the complex phase picked up by the atom travelling 
from u — Z{C, — r) to z = Z{(). This is multiplied with 
the shifted original envelope amplitude x(^t, C ~ ''")• The 
factor U{z,u; k) corresponds to the evolution of the atom 
at rest with the complex mass ^^^jz) (63). From (62) we 
see that U contains, in general, complex dynamical and 
geometric phases. In this way we make the connection to 
[11,12] where such phases were studied for atoms at rest. 

A typical ansatz for the initial envelope amplitude x 
in (65)-(67) is a Gaussian 



x{xT,C)=^f{^^) 

X exp 
ctt, ctl > . 



-3/4. 



-1/2 



4^2^^ 4a2^' 



(68) 



Here A/" > is a normalisation factor, the transverse 
and <7l the longitudinal width, respectively. Supposing 
that at t = we have exactly one atom in this wave we 
get M from 



d^x \i^{x,t;k)\^' 



dV / dz\x{xT,Uz))\' = 1 . 

(69) 



In Appendix B we study the first correction term in 
the series (47) for the function (68). We find that for the 
interferometer parameters specified in (19) ff. the zero or- 
der solution (65) is valid with a 1% accuracy in the region 
of interest provided 



(Tl > 100 /im 
<tt > 100 /im 



(70) 



see (B.64)-(B.66). Under typical experimental conditions 
[1] the actual longitudinal widths of the wave packets used 
are of order 0.4 nm, that is, much smaller than those in 
(70), while ar can indeed be of order 100 /im. Such wave 
packets are described as superpositions of the solutions 
(65) which serve as wavelets in the analysis, see Section 4. 



4 Interference effects 

In this section we shall study, as a concrete example, in- 
terference effects of two 2S states travelling through the 
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interferometer specified in (19) ff. Of course, the formu- 
las given below are easily adapted to other experimental 
setups and more than two states. 

Let us consider the case that we have at time t = 
two narrow wave packets around z « in the interferom- 
eter of (19) ff. These wave packets are supposed to form 
a coherent superposition of two components |a(2:o)) a-nd 
|/3(zo)) of the 2S states where 



a,/3e / C {9,10,11,12} 



(71) 



Here / is defined as an index set with two elements. The 
two waves travel through the interferometer and are sup- 
posed to be analysed &i z = Za as discussed below. For 
the initial state we make the following ansatz 



with 



and 



ael 



(73) 



dk 



^ix,0) = j—f{k)cM^k{z~zo)]x{xT,z~zo). (74) 

Here and in the following all integrals run from — oo to 
-|-oo if not indicated otherwise. This initial wave function 
is supposed to be non zero only in the field free region 
around zq. Taking into account (70) we shall now suppose 
that the parameter SI specifying the field free region in 
(10) has the value 



61 = 10 mm 



(75) 



Note that for zero external fields the initial states |a(zo)) 
with a G I are orthonormal in the ordinary sense since 
these states have different values (F, F3) which are "good" 
quantum numbers for £ = 13 = 0, see Table 2 of Appendix 
A. In (74) x{xt,z — zq) is an envelope amplitude with 
widths in transverse and longitudinal directions 

= 10 mm » (TT, CTL > 100 ptm . (76) 

We suppose the normalisation condition 



^d^ajT Jdz\x{xT,z 



1 



(77) 



to hold, see (69) with 3(z) — z — zq. The Fourier transform 
of X is denoted by 



We require the state in (72) to be normalised which implies 
d^x\<F{x,0)\^ = 



X j d^x'r J dz'cxp[i{k-k')z'] 

X \x{x't,z')\' 
= 1 . (80) 

We define the Fourier transform of !f'(a;, 0) in (74) as 
!^(q'T,gL) 

= J '^^^^T J dz cxp{-iqTXT ~ iqLiz - zo))^{x,0) 
dk - 

^fCkmqT,qL~~k) ■ (81) 

Equations (71) to (81) define the general type of wave 
functions for which we shall calculate interference effects 
below. For later use, we define Gaussian functions as con- 
crete examples for such wave functions: 

x{xT,z-zo)^{27T)-y^a^'a-'/^ 

X exp 



T 

1 . 1 



1/2/^2 ,2x1/4 



fik)^iAnaky^\ai+ai) 

X exp [-crfc(fc - kmf] 



(82) 



(83) 



with 10mm ^ ax, cl > 100 /zm and ^ a/,. We get 
then 



(84) 



xiQT,qL) = (2^)3/42,7^72^ 

X exp(-cr|q| - alql) , 

#(qT,9L) = (2^)3/42aTyK 

X exp[-a^ql - a'^{qL - k„i)^] (85) 

where 



(86) 



Note that for 0"^ <C ctl we have a'j. ~ Ok- 

The evolution of the general state (72) with time is 
obtained from the results of Section 3. Here, of course, 
we have to reinstate the index a on all quantities. Thus 
we have potentials Va, variables Cqi final solutions (65) 
tlja{x,t; k), etc. We get in this way from the initial state 
(72) the solution 



T^fik) Ca'ipaix, t: k)\a(z)) 



(87) 



X x{x't, z') . (78) where 

We superpose functions x in (74) using a function /(fc) 
which is supposed to have a sharp peak at fc = fc„i (w for 
maximum) and a width Ak = with 



ipaix, t; k) ~ exp 

X Ua{z, u\ k) 



i—t + ik{z - zo) 
Zm 



< (Tfc < CTl , kr,i > cr^ 



(79) 



«=2c.(3„(z)-(fc/m)t) 

X xixT,ia{z) - {k/m)t) . 
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Now we come to the analysis of the state (87) at z ~ Za 
where again we have a field free region for z several 
around Za, see (10) and (76). We have, therefore, 



\a{z)) = \a{za)) = \a{zo)) 



(89) 



for \z — Za\ < SI = 10 mm. We suppose that the state (87) 
is projected onto 



\p) = 



where 



This gives 



with 



Wp{x,t)) = \p){p\^{x,t))=%{x,t)\p) 



Y^{a{zo)\pl 



(90) 

(91) 

(92) 
(93) 



%{x,t)^{p\^{x,t)) 

= / '^f{k)Y.pWMx.t\k) . (94) 



We suppose that the total integrated flux Tp in the 
above projection is measured at z = Zq. This is given by 

Tp^ J dt J d^XTjz{xT,Za,t) (95) 

where the z-component of the probability current is 
1 d 

j,{xT,z,t) = —W;(x,t)—Wp(x,t)+c.c. . (96) 



dz 



In evaluating jj, and Tp we take into account (79). We get 
then, as explained in detail in Appendix C, the following 
results 



jz{xT,Za,t) = —'I'pixT, Za,t)>I'p{xT, Za,t) , (97) 



^ a, pel 

X exp[-i(fcs - fc™)(Z\r/3 - Atc)] 

xU*j{Za,Zo;Kn)Ua{Za,Zo]km)^. (98) 

In (97), (98) and in the following only the leading terms 
for large km are kept. In analogy to (56) we define, with 
k = km, the shift of the reduced arrival times at z = Za to 
be 

m f^" 

ATa = 3<.{Za) - (Za " ^o) = TF / '^"(^^ (9^) 



with 0{k„^^) terms neglected, see (C.5) and (C.6). Fur- 
thermore, we introduce the function 



ikd (z' + ^At0 + ^Ato 
X x*ixT,z' + Atp)x{xt,z' + At^) . (100) 



g{ks,ATp,ATa) : 

X / dz' / d^a;Texp 



From (62) we get with k ~ km and fc(zo) = k(za) ~ km, 

Ua{Za, Zo] km) = exp[-i(/)dy„^Q + i(?!>gcom,a] (101) 

where the complex dynamical and geometrical phases are 



777 / ^ 1 

r- dz'[y„(z')--r„(z')] , (102) 

J zq ^ 



Pgcom.Q — Taa {Za)- (^0) 



(103) 



In (102) terms of 0(fc„^) are neglected. From (99) and 
(102) we find 



Re( 



(104) 



In (98)-(104) we have given our general result for the 
integrated flux Tp of the atom in the state |p) where we 
supposed large km, see Appendix C. That is, (98)-(104) 
are valid for all wave packets of the type described in 
(72)~(81). We insert now the Gaussian distributions for % 
and /, see (82)-(86). After a simple calculation presented 
in Appendix C we obtain for this case for Tp, an easily 
accessible observable for lABSE experiments, the result: 

•^P = E exph(Z\T^ - AT^fl{%G'^)\ 

-2^ai ■2'Oi "-rn )\Ja{Za,ZQ;km) ■ (105) 

As a concrete apphcation we discuss our result (105) 
for the following field configuration. We choose the electric 
field as £{z) = and the magnetic field similar to (57) 

B(z)-6(z)e3 , 

B{z) = 10 (e-6(--i-25)^ _ s ■ e-6(--3.75)^) , (106) 

where a parameter s is inserted in order to detune the 
second half of the magnetic field configuration around the 
symmetric arrangement, see Figure 6. The reason for the 
introduction of s will be discussed below. Note that such 
magnetic fields can be produced by suitably shaped coils. 
The so-called spin echo coil produces the scalable part of 
B{z) via controlling the electric current through the coil. 
Figure 7 illustrates the global behaviour of the two states 
a = 9 and a = 11 for different values of s. Setting s = 1 
leads to the roadmap plot shown in Figure 5. Considering 
s = 0.9, the state a = 9 accumulates a delay compared to 
a free wave, while the state a = 11 picks up an advance. 
The opposite is true for s = 1.2. 

Eventually, the evaluation of the total integrated flux 
J-p is shown in Figure 8. In our calculation we choose for 
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Fig. 6. The magnetic field as used for the calculations in this 
section is shown for several values of the parameter s from 
(106). 




Fig. 7. Two examples of the roadmap plots for the states q = 9 
and a = 11 based on the magnetic field (106) and £{z) = 
with fixed parameters a = 0.9 and s — 1.2 respectively. 



the initial wave packet Gaussian momentum distributions 
for the functions in (82)-(85) with 



/c„ = 5.6 X 10^" m-^ 
o'L = o'T = 100 ^m , 
a'l. = 0.4 nm . 



(107) 
(108) 
(109) 



Furthermore, we choose the coefficients Cq of the initial 
state the same as those of the analysing state, namely 



cg = cii = P9 = pii 



1 

71 



(110) 



The initial state is then according to (72) and (81) 



Mx,m,^, = 'Fix,o)-^ (i9) + iii)) 



=0 ' (111) 



d^grdgL 



cxp[iqTXT + iqL{z - zq)] !^(9t, qh) 



(27r)3 

(27r)-3/4aj;i(a^)-i/2 exp[ifc,„(z ~ zo)] 



X exp 



- z 1 / N 



4ct: 



4-^ 



(112) 



Thus, with the choice (110) we superpose the states a = 9 
and a = 11 at zo = and project the evolved state at 
Za = h m onto the same superposition. 




0.7 O.J 



0.9 



1.1 



1.2 1.3 



Fig. 8. The total integrated flux for the magnetic field in (106) 
with the parameter s varied over the interval [0.7, 1.3]. Every 
value of s leads to a fixed flux of analysing states. 



The signal to be observed experimentally is the inte- 
grated flux as function of the parameter s at a fixed point 
on the z-axis, in our case at the final point z^ = 5 m of 
the interferometer, see Figure 8. A variation of s alters the 
terms [Arp — Ato) in (105) as well as the dynamical phases 
0dyn (102). Thus, with the simple field configuration of 
Figure 6 it should be possible to create an interference 
signal simply by tuning the current of the spin echo coil 
and counting outgoing atoms. This signal could serve as 
a reference for more complicated field configurations. 

Let us now consider a configuration with magnetic and 
electric fields as a second example. We define the functions 



6(z) = 10(e-i"(^-°-"5)^+e-i°(^-i-^5) 
— e 



10(2-3.25)'' _ g-10(2-4.25)^ 



(113) 



and 



E{z) = e-30"(-0-9)' + e-30"(-i-6)^ (114) 



where z is given in meters. Now we choose the magnetic 
field B{z), given in /xT, and the electric field £{z), given 
in V/cm, as follows 

fii = 62 = , 

B{z) , z e [0, 0.75] U [1.75,2.5] , 

_ , S(0.75) , z e [0.75,1.75] , (115) 

^■^ ~ s • Biz) , z e [2.5, 3.25] U [4.25, 5] , 

s -6(3.25) , ze [3.25,4.25], 
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and 



£i 



£{z) 
£(0.9) 



z e [0,0.9] U [1.6,5] , 

z £[0.9,1.6], (116) 



f 2 = £^ = 



respectively. In Figure 9 we show B^{z) and £i{z). The 
electric field is ramped to a constant value and back again 
only within a region of constant magnetic field. This pro- 
cedure implies vanishing geometric phases. 



H 
a. 




> 



z [ni] 

Fig. 9. The magnetic field component Bs{z) (115), solid line, 
and the electric field component £i{z) (116), given as the 
dashed line. 



For the calculation of the spin echo signal with the 
external fields (115) and (116), we again employ the initial 
specifications (107)-(112). The results for both the pure 
magnetic field and the combined magnetic and electric 
fields arc shown in Figure 10. In the presence of the electric 
field 5, we find a decreased total integrated flux J^^^°(s) 
due to the enhanced decay rate, compared to the case of 

In order to interpret Figure 10, we now give an es- 
timate of the quantities which lead to the altered spin 
echo signal when the electric field (116) is switched on. 
For 5 7^ the quadratic Stark effect leads to a shift of 
the complex energies Ea. The shift of the real parts Va 
of Ea turns out to be negligible. Therefore, the reduced 
arrival times do not change substantially compared to the 
case of 5 = 0. The change of the oscillation frequency 
of J^p^°{s), induced by Va, is also not significant if we 
compare to J-^^^{s), see Figure 10. However, the altered 
imaginary parts — have a large impact on the total in- 
tegrated fiux via the dynamical phases. For the case of the 
non-zero electric field (116), wc get for the states a = 9, 
10 and 11 the relation 



rf'^°(z)«30ori^=°)(z) 

« 300 Ts w 2470 s-i 



(117) 



for the parameter range 0.85 z < 1.65, where Eg is 
given in (23). Taking into account that in this case the 
geometric phases are zero and considering (105) together 



0.8 



0.6 



0.4 



0.2 




0.8 



0.9 



1.1 



1.2 



Fig. 10. The spin echo signal corresponding to the case where 
only the magnetic field (115) is present is given by the (red) 
dashed line. The spin echo signal corresponding to the presence 
of the magnetic field (115) and the electric field (116) is given 
by the solid line. 



with (101) and (102), we find for s = 1: 



^^^^o)(s = l)«cxpi 



■300 



2 — 1.65 m 



2^0.85 m 



dz'Es) 



Inserting 



we get 



(3500m/s)~i , 



(118) 



(119) 



Tf^°\s = 1) « 0.569 J'^^^=°)(s = 1) « 0.562 . (120) 

The rough estimate (120) for the total integrated flux co- 
incides satisfactorily with the numerical result 



1) w 0.569 



(121) 



given as the global maximum of the spin echo signal for 
£ 7^ in Figure 10. With this we close our discussion of 
the concrete examples. 

In future work we will study field configurations where 
geometric phases, both P-conserving and P-violating ones, 
occur. In principle the spin echo experiments should pro- 
vide the experimental access to these geometric phases 
since a non-vanishing geometric phase will modify the sig- 
nals like that of Figures 8 and 10 for which only the dy- 
namical phases are non-zero. For example, switching on 
suitable electric fields, geometric phases will contribute 
to the total phase factors f/| and Ua in (105). Besides 
the changed magnitudes of the peaks, the crucial infor- 
mation is the displacement of their positions in the s — 
J'p{s) interference-diagram compared to the reference sig- 
nal indicating non-zero geometric phases. Of course, a P- 
violating Berry phase will modify the reference signal in 
a slightly different way if the measurement is performed 
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with the space reflected field configuration. In this way, the 
high sensitivity of lABSE experiments could possibly lead 
to measurements of P-conserving as well as P-violating 
geometric phases in mctastable hydrogen. 

5 Conclusions and outlook 

In this article we have discussed the evolution of the mcta- 
stable 2S states of atomic hydrogen within the setup of 
longitudinal atomic beam spin echo (lABSE) experiments. 
The conditions for adiabatic evolution and sufficient flux 
at the end of the interferometer were investigated. Using 
the results for atoms at rest [12], we derived the solution 
of the effective Schrodinger equation for the mctastable 
atoms for the case that the longitudinal-kinetic energy 
is much larger than the transverse-kinetic and the po- 
tential energies. The wave function consists of a plane 
wave factor, a complex phase factor containing a geomet- 
ric Berry phase and the envelope amplitude of the prop- 
agating state. In Section 4 we considered the interference 
effects between two mctastable 2S states and computed 
the total integrated flux J^p for a given projection at the 
end of the interferometer in (105). This quantity J-p is 
conveniently accessible to experiments and encodes the 
information from the phase differences picked up by two 
mctastable 2S states in the interferometer. In future work 
the extraction of geometric phases from such experiments 
will be further analysed. The main result of the present 
paper is to provide the suitable framework for the rigorous 
theoretical description of stable or mctastable atoms in a 
longitudinal spin echo device. Here we have considered 
the adiabatic case and abelian geometric phases. Possi- 
ble extensions of the present work would be to consider 
non-abelian dynamic and geometric phases as well as non 
adiabatic motion. 
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Appendix 

A Values for quantities related to n = 2 
hydrogen 

In this appendix we collect the numerical values for the 
quantities entering our calculations for the hydrogen states 
with principal quantum number n = 2. We specify our 
numbering scheme for these states and give the expres- 
sions for the mass matrix at zero external fields, and for 
the electric and the magnetic dipole operators. 



In Table 1 we present the numerical values for the 
weak charges x = 1, 2, the quantities Aq, the Lamb 

shift L — i?(2Si/2) — E{2Pi/2), and the fine structure 
splitting A = £'(2P3/2) - E{2Pi/2)- Here Aq denotes 
the total polarisation of the nucleus carried by the quark 
species q {q = w, d, s). The ground state hypcrfinc split- 
ting energy is denoted by A. We have A = E{lSi/2,F = 
1) — E{lSi/2,F = 0) for hydrogen. For the precise defini- 
tions of all these quantities see [12]. The n = 2 states of hy- 
drogen in the absence of P-violation and for zero external 
fields are denoted by \2Lj,F,F3), where L, J, F and F3 
are the quantum numbers of the electron's orbital angular 
momentum, its total angular momentum, the total atomic 
angular momentum and its third component, respectively. 
The quantum numbers S for the electron spin and / for 
the nuclear spin are omitted, since these are fixed quan- 
tities. In the following the ordering of the atomic states 
in the matrix representations of operators is according to 
Tabic 2 where we give the numbering scheme for the states 
which we consider. For electric field £ and magnetic field 
B equal to zero we have the free 2S and 2P states. We 
write L, P, S since these states include the parity mix- 
ing due to Hpy, see (1) of [12]. Consider first atoms in a 
constant B-field pointing in positive 3-direction, 

B^Seg, B>0. (A.l) 

The corresponding states \2Lj, F, F3, 0, Be^) are obtained 
from those at S = by continuously turning on B in the 
form (A.l). Of course, for \B\ ^ 0, F is no longer a good 
quantum number. Here it is merely a label for the states. 
We now choose a reference field Srcf = Brcfe-i, B^of > 0, 
below the first crossings in the Breit-Rabi diagrams, for 
instance B^cf ~ 0.05 mT. We define the states |2Lj, F, F3, 
£, B) for arbitrary B fields in the neighbourhood of B^cf as 







Ref. 


z 


1 




N 







I 


1 
2 






0.04532(64) 


(11) of [12] 


Si{Z,N) 


-2.78(4) • 10"^^ 


(20) of [12] 


Au{Z,N) - Ad{Z,N) 


1.2695(29) 


[18] 


As{Z, N) 


0.006(29) (7) 


[19] 


q'w\z,n) 


-0.1145(31) 


(12) of [12] 


&2{Z,N) 


7.04(19) • 10"^^ 


(20) of [12] 


L{Z, N)/h 


1057.8440(24) MHz 


[20] 


A{Z,N)/h 


10969.0416(48) MHz 


[20] 


A{Z, N)/h 


1420.405751768(1) MHz 


[21] 



Table 1. Values of parameters for numerical calculations. 
The weak mixing angle in the low energy limit, sit? dw ~ 
0.23867(16), was taken from [22]. The uncertainty in 5i is dom- 

2 (2) 
inated by the uncertainty of sin -dw- The uncertainties in 

and ^2 for hydrogen are resulting from the errors of the weak 

mixing angle and As in equal shares. 
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the states obtained continuously from |2I/,/, 0, Brcf ) 

by turning on £ and [B — Brof) as \£ and \{3 — Sicf), 
respectively, with A G [0, 1]. Then both F and are only 
labels for these states and no longer good quantum num- 
bers. 



hydrogen 



a \2Lj,F,F3,S,B) 









1 I2P3/2 


2, 




2 I2P3/2 


2, 


i,e,B) 


3 I2P3/2 


2, 


0,£,B) 


4 I2P3/2 


2, 


-l,£,B) 


5 I2P3/2 


2, 


-2,£,B) 


6 I2P3/2 


1, 


1,£,B) 


7 I2P3/2 


1, 


0,£,B) 


8 I2P3/2 


1, 


-l,£,B) 


9 1 25*1/2 


1, 


l,£,B) 


10 |2S'i/2 


1, 


0,£,B) 


11 |2Si/2 


1, 


-1,£,B) 


12 |2S'i/2 


0, 


0,£,B) 


13 I2A/2 


1, 


1,£,B) 


14 I2A/2 


1, 


0,£,B) 


15 I2A/2 


1, 


-1,£,B) 


16 I2A/2 


0, 


0,£,B) 



Table 2. The numbering scheme for the atomic states of hy- 
drogen. 



Tables 3-5 in Appendix D show the non-zero parts of 
the mass matrix ^p((5i , S?) of (4) for hydrogen for zero 
external fields, of the electric dipole operator D and of 
the magnetic dipole operator /x for the n = 2 states. 

Next we discuss the conditions for decoupling of the 
evolution of the 2S states from the 2P states. We use for 
that the results of Appendix B of [12]. Let 

T « 1 ms (A.2) 

be the total observation time. We then have, with Vz from 
(20), v^T « 3.5 m which is of the same size as I (19). From 
(B.1)-(B.3) of [12] we get as conditions for the decoupling 

2C12 



TAF^lr, 
4Cl2C21 

2C21 



< 1 



< 1 



< 1 



TAP 
with 

^-fniin ~ Fp = Tp^ . 

In the following we shall use with rp from (23) 
VzT w Z w 5 m , 
VzTp = 5.6 /im . 



(A.3) 
(A.4) 
(A.5) 

(A.6) 
(A.7) 



As in [12] we shall use for our estimates the quantities £q 
and Bq which characterise the values of the field strengths 
where the electric and magnetic energy, respectively, reach 
the value of the Lamb shift. With ^3(1) the Bohr radius 
of the H atom and fj,B the Bohr magneton we have (see 
(B.IO), (B.19) and (B.20) of [12] 



£0 



L 



477.3 V/cm , 



= 43.65 mT 



(A.8) 
(A.9) 



For a time- varying electric field £{t) the quantities C12 and 
C21 are estimated in (B.5)-(B.9) of [12]. We take this over 
for our case with the replacement dt v~^dz, that is. 



d£{t) 



dt 



d£(z) 



In this way we get 



C12 « C21 



Tvz 
V3£o 



■ max 



dz 



d£{z) 



dz 



(A.IO) 



(ATI) 



Here and in the following the maximum is to be taken for 

zo < z < Za- We set 



max 



d£{z) 



dz 



AE_ 
~A'z 



with a typical variation (see (25)) 
A£ w IV/cm 



(A.12) 



(A.13) 



and Az the allowed variation length which we want to 
estimate. Inserting (A.12) in (A. 11) we get from (A.3) 
and (A.5) 

2 A£ 

Z\z > —w^Tp— w 0.02 /^m . (A. 14) 

V3 to 



From (A.4) we get 



-AS 



Az > -^^JvzTvzTp— w 12/im 
V3 '^o 



(A.15) 



Thus, for electric field variations (A.12), (A.13) a length 
scale 



Az > 50 /im 



(A.16) 



should be sufficient to guarantee decoupling of the 2S from 
the 2P states. 

The allowed time variation of magnetic fields was es- 
timated in (B.17)-(B.22) of [12] but allowing rather large 
electric fields of order £q/2. Here we allow only fields of 
order 



£ = IV/cm 



(A.17) 
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Thus we get from (B.5) and (B.18) of [12] with the inser- 
tion of an additional factor 2£ /Eo and the replacement 



dB{t) 



dt 



d3{z) 



dz 



the estimate 



2 E 

ci2 ~ C21 ~ —;=Tvz max 



We set here 



a/3 £oBo 
dB{z) 



dB{z) 



dz 



(A.18) 



(A.19) 



dz 



AB 



Z\B = 1 mT . 

Inserting now (A.19) in (A. 3) and (A. 5) we get 

4 £AB 
V3 fo^o 

From (A. 4) we get 

4 



6 X 10 ^m 



Az > -^\/ v^Tv^Tp^^- K. 0.6 /xm 
V3 cooo 



(A.20) 



(A.21) 



(A.22) 



Thus, a variation length Az w 10 /xm should be sufficient 
to ensure decoupling of the 2S from the 2P states for the 
magnetic fields considered. 



B Solution for one component, details 

In this appendix we give the detailed derivation of the re- 
sults of Section 3. First we derive the invariance of Va{z) 
in (32) under the local phase transformations (34) where 
the functions r]a{z) are assumed to be real and differen- 
tiablc. From (29) we obtain 



d d 

^Jaa{z) = i{a{z)\ — \a{z)) . 
oz oz 

Hence the transformation (34) yields 



(B.l) 



(B.2) 



and 



d d 

Ko'iz)\j-\a{z)) - j-Va{z) 
oz oz 



OZ 



-2i(^^77„(z)^ (a(z)l^la(z)) 



(B.3) 



This implies invariance of the potential Va{z): 
Va^iz) V'^iz) ^ E^iz) + Eo 



1 

2m 



d 

{a{z)\ — \a{z)) 



2i ( §;^o.iz)j ia{z)\-^\aiz)) + (-^V^i^) 



d \ ----- Q 
2i( ^^?a(^)j {a{z)\ — \a{z)) 



d_ 

^dz 

Vaiz) 



(a(z)l 



d_ 

dz 



a{z)) 



(B.4) 



Next we discuss the derivation of (45), the equation for 
the envelope function A, and of the solution (47). We use 
the variables r, C defined in (38) and make the ansatz (40) 
for the wave function with, in general, complex 0(z,t). 
Inserting (40) in (36) gives 



2i m 



dA 90^4 
dt dz dz 



2m 



dcj)^ 
~dt 



d(j)_ 
dz 



2mV{z) + im_r(z) 
.d\ 



dz'^ 



A 



(B.5) 



The next step is to insert </> from (43) and to change vari- 
ables from XT,z,t to xttC,,t, see (38), (39). Some useful 
relations in this context, with k{z) defined in (42), are 



d_ 
di 

d_ 

dz k{z) d( 



k d 
m dr ' 
k d 



dk{z) _ TO dV{z) 
dz k{z) dz 

Furthermore we have 

dr = —dt , 

TO 



(B.6) 



(B.7) 



dC 



Hz. 



-dz 



(B. 



The local longitudinal velocity of the atom at position z 
is 



v{z) 



k{z) 



(B.9) 



Suppose now that the atom arrives at position z at 
time t. Then it has to start at a certain position u{z, t) at 
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time t = 0. Wc have 



dz = w(z)dt, 



dz' 

u(z.t) V{Z') 

u 

dz'- 



dt' = t 



■,(z,t) Hz') m 
We can rewrite this as 



t = T 



u{z^t) 



dz' 



k{z') 



dz'- ^ 



fc(z') 



(B.IO) 
(B.ll) 

(B.12) 
(B.13) 



Thus u{z,t) is a function of C — t only. We can express 
this with the help of the function Z(-) defined in (39): 



u{z,t) ^ Z{C - t) ^ ZOiz) - -t) 

m 



(B.14) 



Expressed in the r, ^ coordinates we have the following. 
The atom arrives at reduced time r at position Q if it starts 
at T = at position C ~ From (B.ll) wc get easily 



du{z,t) k{u{z,t)) 

dt m 
du{z,t) k{u{z,t)) 



dz 



k{z) 



(B.15) 



Using this we find for the partial derivatives of (f> (43) the 
following: 



d<t>{z,t) 
dt 



fc2 i 1 dViu) i , , 



Here and in the sequel k{u) is to be read as k{u{z, t)), r{u) 
as r{u{z, t)), etc. 



d<f>iz,t) im ( 1 dV{z) 



d^(f>iz,t) 



2k{z){k{z) dz 
-{z^u) 

m dV{z) 



(B.17) 



k{z) dz 
m? dV{z) 
'2fc3(z) dz 



1 dV{z) 
k{z) dz 

— {z —f u) 



r{z) 



2fc2(z) 

d^Viz) 
dz-^ 



m fdV{z) 
kHz) V dz 



kiz)^^-{z^u) 



(B.18) 



Inserting (B.15)-(B.18) in (B.5) we get the following 
equation for A{xt, Ci ''')■ 



^§i + Hz)-^]AixTX,T) 



1 J 

2 j4fc2(z) 

dV{z) 
2fc3(z) dz 



1 dV{z) 



k{z) dz 



— r{z) — {z —t u) 



1 dV{z) 



k{z) dz 



-r{z)-{z^u) 



m 



2k^{z) 



dv{z)Y d^v{z) 



k'^{z) \ dz J dz'^ 
A{xt,C„t) 



2k{z) 



1 dV(z) 
k{z) dz 



— r(z) — (z u) 



dA. . . 
^i-T,Cr) 

(B.19) 



Replacing the derivatives d/dt and d/dz by d/dr and 
d/dC according to (B.6) leads to 



d d 



AixT,CT)^{LA)ixT,C,r) . (B.20) 



dr dC 

Here L is the following differential operator 

d d'^ i 

L = loiz,u) + hiz,u)— + hiz)-^ + ^At , (B.21) 



dV{z) 



2fc3(z) dz 



1 dVjz) 
k{z) dz 

1 dVjz) 
k{z) dz 



r{z)-{z^u) 
-r{z)-{z^u) 



m f dV{z) V d^Vjz) 
fc2(z) \ dz J dz"^ 



- k{z)^^ - {z ^ u) 



dz 



h{z,u) 



2fc2(z) 



1 dVjz) 
k{z) dz 



r{z) ~{z-^u) 



Hz)^i 



m dV{z) 
2fc3(z) dz 
k 



2fc2(z) 



(B.22) 

(B.23) 
(B.24) 



In (B.21)-(B.24) we have to consider z and u as functions 
of C and T according to 



z = Z{C) , 
u^ZiC^r) 



(B.25) 
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see (39) and (B.14). We arc interested in the case where where we define the operator K as follows: For any func- 
(6) holds which implies also tion /{xt, Ci ''") we set 



fc^ > |2mV(z)| . (B.26) ^ n /"' A^'frfsr 'w 

- ^ C^"^ 

Then k{z) w k and we find that all terms of L have at 3g^ 

least one power of 1/k. Thus, for large k we get by formal 

power comiting Since L (x 1/k, see (B.27), we have also 

1 K(xl/k. (B.37) 

L cx ^ . (B.27) 

That is, all terms of K have an explicit factor 1/k. The 
The next step is to define the Green's function for solution of (B.35) is easily obtained 
(B.20): 

^=(l-if)-V = E^>- 

— + \GiCr)=SiOSir). (B.28) 

We write this as 

One solution is _^ ^ ^^(o) ^ ^^(i) _^ ^ _ ^ 39) 

G(C, r) = e{C + t)S{C - t) . (B.29) where 

With (B.29) we can transform (B.20) into an integral A''^\xt,Ci''') ~ ^fii^TX ^ t) , 

Aixr, C, r) = Bixr, C - r) + f dC f dr' ^'"^ = C, r) , 

xG(C-C',r-r')(L^)(=rT,C',T') • • (^-40) 

= B{xt,C — t) The successive terms in (B.39) have higher and higher 

/(^ powers of according to (B.37) and are thus more and 

d(' {LA){xT,C\'''')\r'-T c+C ■ more suppressed for large k. This is, of course, only a 

'oo formal argument. With (B.39) and (B.40) we have given 

(B.30) the proof of (45)-(47). The zero-order term is further dis- 
cussed in Section 3. The range of validity of the expansion 

Here B{xt, C ^ ''') satisfies the homogeneous equation (j^ 39") jg discussed below 

/ ^ „ s Let us consider now the first-order term in (B.40) 

— - + — I B(xT C-t)^0 (B 31) for a general function ip. With (B.21)-(B.25) and (B.36) 

dCj we get 

and is determined by the initial condition at r = 0. Indeed, 
let us define 



0'^ 

^ixT,0---AixT,Cr)Uo ■ (B.32) + hm'), ZiC - r'))— + UZ{C))g^ 



We set 



B{xtX-t) = ip{xTX-T) fw^/^/x ^/^ , , /^.WN ^^ 9 



'J —00 



= |z„(Z(C'), Z{C - r)) + /i(Z(C'), ZiC - r))^ 

MZiC))^ + ^^t| ^{xt, C-r), (B.41) 



and obtain from (B.30) 

AixTX^-r) = ip{xTX-T) A^^\xtX,t) = iKA'^°^)ixTX,T) 



f dC'(L^) (a;T,C',T')l.,.._c+C 



^ dC{LA^"^){xT,C,r') 
(B.34) r d 

Ko{C,t) + «^i(C,t)^ + «^2(C,T")^2 



We write this in matrix notation as 

A = ip + KA, (B.35) ' 2k 



AT^A^'\xTX,r) . (B.42) 
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Here we define 
«^o(C,t) 

'«i(C,t) 



c 



/ dc';i(z(c'),2(c-r)) 

JC~T 



(B.43) 



C-r 



Now we shall make rough estimates of the size of the func- 
tions |kj(C,t)|, (j — 0, 1,2). For the parameter r to con- 
sider we suppose as in (A. 7) 



T « Z = 5 m . 
For the potential we take as an estimate 
\V[z)\<V^ , 

= lOV/cm , 
= 1 mT . 

With 

fc = = 5.6 X 10^° m"^ , 
see (20), this gives 

\Viz)\<Vra^lO~^eY, 



\V{z)\ 



2m 



=: 7/ = 1.5 X 10"^ . 
For the variation of the potential we estimate 



dV{z) 



dz 



,^.Ae AB 



AS = 1 V/cm 
AB = 0.1mT , 
Az = 100 ^m . 



This gives 



dViz) 



dz 



< 



lOAz 



Similarly we estimate 

d^Viz) 



dz'^ 



< 



Vrr, 



lOAz^ 



(B.44) 



(B.45) 
(B.46) 



(B.47) 



(B.48) 



(B.49) 



(B.50) 



For the decay rate we estimate according to (24) and 
(B.44) 



r{z) 
drjz) 

dz 



< 



< 



I ml 
k 

mlAz 



(B.51) 



Now we can estimate Kj{(,T), [j = 0,1,2), from (B.43) 
and (B.22)-(B.24) as follows 



I«o(C,t)| 



< 



2k 4fc2 



2v;„ 



2fc3 IQAz 



2y„: 



m 



lOkAz 

2 



\QkAz 

2k 
ml 

V„ 



2k 
ml 



m / 

¥ \10Az J ' lOAz 



2k [A [lO Az 
' 10- , 

I ( m 



3 ?/ 



mlAz 
V 



2 lOlAz 20Az^ 



1 



2kAz 
\^i{Cr)\< 



1 

jAi 
(B.52) 



MCr)|<^ 



k y k 
i_ n 
k 
I 



1 K, 



k IQAz 
3 V 



k 

ml 



1 m Vm 



2 fc2 IOA2 



I 4 10Z\z 



(B.53) 
(B.54) 



With (B.52)-(B.54) we can estimate the size of A'^^^ rela- 
tive to from (B.42). We choose as a concrete example 
for the amplitude the one corresponding to the Gaus- 
sian function x in (68). From (66), (68) and (B.40) we get 



A^^'^XtX^t) ^v{xt,C-t) =AA(2^)-3/V^ia 



-1/2 



X exp 
X exp 



, 2mV{z') . ■ 
dz ^— : — i7(u) 



k{z') + k 



u=Z{C-r) 



Here we suppose 

100 /im < oltOt < 300 /im . 



(B.55) 



(B.56) 



The Gaussian function in (B.55) is nonzero practically 
only for 

\C-t\<2.(7l. (B.57) 

Then we get also for w = Z(C — r) from (38) and (39) 

|w- zol < ScTL ■ (B.58) 

This means that u in (B.55) varies only in the field-free 
region around zq- Thus, the last exponential function on 
the r.h.s. of (B.55) can be replaced by 1 and, similarly, 
we get N = 1 from (69). The amplitude A^'^^ and its 
derivatives with respect to C read then 



X exp 



T "L 
1 2 1 

'4^^^ 4^ 



(C 



(B.59) 
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dC 



'xt,C,t) 



2al 



2al 



t) , (B.60) 



(C 



2al 



(B.61) 



want. The actual narrow wave packets are formed as su- 
perpositions of the wavelets as discussed in Section 4 and 
Appendix C. Finally we emphasise again that we have 
given here only an example of accuracy estimates for the 
wavelet solutions of (B.20). For other experimental condi- 
tions, that is other electric and magnetic fields, other atom 
velocities etc., similar estimates can easily be performed. 



T \ T / 



(B.62) 



Inserting (B.59)-(B.62) in (B.42) we get 



a^^\xt,(:,t) 



ko(C,t) 

'«2(C,t) 

ir 1 
'2k2^ 



2al 



iC-r) 



2 1 



2,| 



^j}.4'".,x.,C.). 

(B.63) 



In the only interesting range given by (B.57) we get with 
(B.52)-(B.54) as an estimate of the order of magnitudes 



A^'H^tX^t) 



< 



MC,r)\ + \K,iC,T)\ 



2ai 



l«2(C,r)| 



2ol 



I 1 
2k a\ 



10" 



2kuL 



2ka\ 



(B.64) 



where the quantity x\l{2a\') from the last term on the 
r.h.s. of (B.63) is neglected due to the exponential decrease 
oi A^^'^{xtX-,t) for increasing \xt\- Numerically the last 
two terms on the r.h.s. of (B.64) are the dominant ones 
and we find 











- 2k \2al 





(B.65) 



Inserting here /, k and (Tl, cjt from (19), (B.46) and (B.56), 
respectively, we get 



2fc 



1 



2al 



< 0.01 



(B.66) 



To summarize: we find that with our conditions speci- 
fied above the zero order amplitude A^'^^ is accurate with 
errors of at most 1%. The origin of the dominant errors 
shown in (B.65) is easily identified. It is due to the fact 
that A^'^'' docs not include the longitudinal and trans- 
verse dispersion of the wave packet. But this presents no 
problem, since we have discussed in Section 3 and in this 
appendix only wavelets, where we can - within reason- 
able limits - choose the widths cfl and ax as large as we 



C Interference effects, details 

Here we represent the details of the calculations of the in- 
terference effects discussed in Section 4. The z component 
of the probability current for the wave function \tljp{x,t)) 
(92) is given in (96). Inserting there (94) and taking into 
account that the k distribution is supposed to be sharply 
peaked at km, see (79), leads to (97). This gives for the 
integrated flux Tp (95) the result 



^rCk')fCk) 



(27r)2 



X exp 



^2 _ ^/2 

-i — 1 + i{k - k'){za - zq) 

Zm 



X Ua{Za,U]k)\ 



k' k 

X X*{xT-,Ci3 t)x{xT-,Ca 1) 

m m 



where 



Now we make the following approximations 



) 

(C.l) 
(C.2) 



1 — - (fc - ~k')^^-^t « (fc - k'Y^t 

2m 2m m 



h 



exp 



with 



2m 



-t + i(k-k'){Za- Zo) 



Zo 1 . 

m 



exp [i{k — k')z''\ 
(C.3) 

(C.4) 



Furthermore we have 



Cq - — t = 3aiZa) - — t = ATa{Za; k) + Z<j 

m m 



k 

—t . 

m 

(C.5) 



Here ATa{za;k) gives according to (56) the difference of 
reduced arrival times of the wave (88) relative to the free 
wave. We have with (38), (39) and (61) 



ATaiZa;k) = dz' 

Jzo \ka{z) 



1 



dz'- 



2mVa{z') 



k^{z'){k^{z') + k) 



(C.6) 
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The function 



X XT,Ca 1 

m 



is significantly non zero only for 



Ca 1 

m 



(C.7) 



(C.8) 



(C.9) 



varies only over a distance of a few ctl around zq, that is, 
in the field-free region 



with (Ti = 100 to 300 fini. Correspondingly, 



U — Za [ Ca t 



We consider now (see (62)) 



(C.IO) 



Uaiza, u; k) = exp 



. , , 2mVa{z') 1 ka{Za) 
1 / dz - — r-: — - In ■ 



1 , ,mra{z') 



dz' 



Ka [Z ) 

(C.ll) 



for u given by (C.9), (C.IO). We have then 

ka{Za) — ka{u) = k , 



, 2mVa{z') 
" ka{z')+k 



ka{u) 



, 2mVa{z') 
' ka{z') + k ' 



laa{u) = laa{zo) 

Also the error made by setting 



dz 



,mra{z') 



dz 



,mra{z') 



ka{z') J ka{z') 

is negligible. Therefore we can set 

Ua{Za,U;k)\^^^^^^^_j^^-^ = Ua{Za,Zo;k) 

, 2mVa{z') 



(C.12) 



(C.13) 



(C.14) 



(C.15) 



exp 



-i / dz' 

' 20 ka \Z) + k 



(C.16) 



Now we remember that /(fc) in the integrand in (C.l) is 
sharply peaked for k — km- We expand therefore the terms 
in the argument of the exponential function in (C.16) 
around km- For the integral over Vq,(z') which gives a fast 
oscillating term we expand up to first order in (k — fc,„). 



For the remaining terms we keep only the zero-order con- 
tributions. We get then 



, 2mVa(z') 



dz' ^ - = / dz 

kc,{z') + k J,^ k^[z') + k 



, 2mVa{z') 



k=k„ 



- [k -km) dz' 



2mVa(z') 



fc„(z')(fc„(z') + fc) 



k=k„ 



0[{k-kmf] 

, 2mVa{z') 



dz 

zo k(^{z') -\- k 



{k km)^Ta 



k=k„ 



+ 0[{k-kmf] 

where we set 



(C.17) 



(C.18) 



see (C.6) and (99). Inserting (C.17) in (C.16) we find with 
the approximations discussed above. 

UaiZa,Zo;k) « Ua{z 

(C.19) 



Ua{za, Zo; km) ~ exp 

1 



,mVaiz') 



, ,mra(z') . , . . / s 

„ , dz ^ h IJaa(Za) - 17aa(^o) 

^ ./^n f^m 



exp[-10dyn,Q + 1</ 



(C.20) 



Here we use the definitions (102) and (103). Putting now 
everything together we find from (C.1)-(C.5) and (C.19), 
(C.20) 

^P^J ^nk')fik) J dz'J d^z^cxp [i(fc-fc')z'] 
X ^ pi3p*^c*pCa exp [ - i(fc' - fc,„)Zir/3 -f- i(fe - km)^Ta] 

X X*{xt,z' + ATf3)x{xT,z' + ATa) 

^0: '^ni ^ai ^Oi f^m ) ■ (C.21) 

Here and in the following all integrals run from — oo to oo. 
Now we change variables and set 

k(^ — k ~\~ — k . 
2 2 ' 



few 



k-k' 



(C.22) 



This gives (98) with g(ks, At/^, At^) as defined in (100). 
Inserting there the Fourier transform of x according to 
(78) we get 

g{h,ATi3,ATa) 

= J dVdgW9Lr(fc.-^('?L-9L)) 

X f{ks + \{q'L - 1l)) exp [i(ql + qL)\{AT^ - Arp)] 



X X*{qT,q'L)xiqT,qL) - 



(C.23) 
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Finally, wc insert the explicit Gaussian functions /(•) and 
x(-) from (83) and (84), respectively, in (C.23) and get for 
this case 



-2(J^(ks ~ fcm) ^ 



X exp 



8-i 



(C.24) 



With (C.24) we get for the integrated flux Tp from (98) 
the final result (105). 



T. Bergmann et al: Parity Violation in Hydrogen and longitudinal Atomic Beam Spin Echo I 



21 



D Matrices for n = 2 states of hydrogen 

Tables 3, 4 and 5 show the non-zero parts of the mass matrix i^Q{5i, 82) for zero external fields, of the electric dipole 
operator and of the magnetic dipole operator fi for the n = 2 states of hydrogen. We give all these matrices in the 
basis of the pure 2S and 2P states, that is the states for zero external fields and without the P-violating mixing. Thus, 
the P-violation parameters 61,2 occur in the matrix ^^((Si, 62) explicitly. 
In Tables 4 and 5 wc use the spherical unit vectors, which arc defined as 

60 = 63, e± = (ei ±i62) , (D.l) 

where 6; {i = 1, 2, 3) arc the Cartesian unit vectors. For e±, the following relation holds: 

el - -6^ . (D.2) 



Table 3. The mass matrix , ^2) (4) for the n — 2 states of hydrogen for the case of zero external fields. For the explanation 

of the variables A, L and A and their numerical values see the introduction of Appendix A and Table 1. The PV parameters 
5i,2 can also be found in Table 1, the decay rates -Tp.s are given in (23). 





2^3/2,2,2 


2^3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


2Si/2,l,l 


2-P3/2, 2, 2 
















2^3/2,2,1 
















2P3/2,1,1 








96 


A 

192 \/2 





2Pi/2,l,l 








A 
192 ^/2 


96 2^ ^ 




2Si/2,l,l 













^ + ^ 



(Table 3a) 
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2^3/2,2,0 


2A/2,1,0 


2Pi/2,l,0 


2Si/2,l,0 


2Pi/2,0,0 


2Si/2,0,0 


2^3/2,2,0 



















2^3/2,1,0 





A - A 


A 
192V^ 











2Pi/2,l,0 





A 
192 v/2 


96 2 P 


i5iL 








2Si/2,l,0 








-i5lL 


+ H 








2 A/2, 0,0 














A i r 

~32 - I'^P 


i(5iL 
+|i52L 


2Si/2,0,0 














-iJiL 
-|i52L 


r 3^1 
-'^ ~ "32" 



(Table 3b) 





2P3/2,2,-l 


2P3/2,1,-1 


2Pi/2,l,-l 


2Si/2,l,-l 


2P3/2,2,-2 


2P3/2,2, 


-1 
















2P3/2,1, 


-1 







A 

192^2 








2Pl/2,l, 


-1 





A 
192^2 


96 2^ P 


iSiL 





2'S'i/2, 1, 


-1 








-iSlL 


^ + ^ 





2P3/2,2, 


-2 














-^rp 



(Table 3c) 
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2P3/2,2,2 


2^3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


251/2,1,1 


2^3/2,2,2 
















2P3/2,2,1 
2P3/2,1,1 

2Pi/2,l,l 
2Si/2,l,l 





3e+ 











-\/|eo 






— \/3eo 


-•v/feo 




(Table 4a) 





2^3/2,2,0 


2^3/2,1,0 


2Pi/2,l,0 


25i/2,l,0 


2Pi/2,0,0 


2Si/2,0,0 


2-P3/2,2, 1 











3 








2^3/2,1,1 
















-^/6e_ 


2Pi/2,l,l 


















251/2,1,1 






















(Table 4b) 












2^3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


25i/2,l,l 






2^3/2,2,0 

















2P3/2,1,0 











/^-^ 






2Pi/2,l,0 











\/3e+ 






2Si/2,l,0 






y3e+ 









2Pi/2,0,0 











-\/3e+ 






2Si/2,0,0 







-\/3e+ 













(Table 4c) 










2^3/2,2,0 


2P3/2,1,0 


2Pi/2,l,0 


25i/2,l,0 


2Pi/2,0,0 


2Si/2,0,0 


2^3/2,2,0 











\/6eo 








2^3/2,1,0 



















2Pi/2,l,0 



















25i/2,l,0 


V6eo 
















2Pi/2,0,0 











— \/3eo 








25i/2,0,0 





\/6eo 














(Table 4d) 
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2^3/2,2,-1 


2^3/2,1,-1 


2Pi/2,l,-l 


25i/2,l,-l 


2P3/2,2,0 











V 


2^3/2,1,0 











-^^^e 


2Pi/2,l,0 













2S'i/2,l,0 






-^e_ 





2Pi/2,0,0 













25i/2,0,0 





V6e_ 


-y/ie- 






(Table 4e) 





2^3/2,2,0 


2P3/2,1,0 


2Pi/2,l,0 




2Si/2,l,0 2Pi/2,0,0 2Si/2,0,0 


2^3/2,2,-1 




















2^3/2,1,-1 















y|e+ 


-\/6e+ 


2Pi/2,l,-l 















y3e+ 


\/3e+ 


2Si/2,l,-l 












^/3e+ 










(Table 4f) 












2P3/2, 2, 


-1 2P3/2,1, 


-1 2Pi/2,l, 


-1 2Si/2,l,-l 


2P3/2,2,-2 




2^3/2,2,- 


1 




















2^3/2,1,- 


1 




















2A/2,1,- 


1 




















2Si/2, 1, — 


1 




y^eo \/3eo 





3e_ 




2^3/2,2,- 


2 













— 3e+ 








(Table 4g) 
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Table 5. The suitably normalised magnetic dipole operator /x//is for the n = 2 states of hydrogen, where ^ib = eh/(2me) is 
the Bohr magneton and g — 2.002319304(76) is the Lande factor of the electron [23]. 





2^3/2,2,2 


2^3/2,2,1 


2^3/2,1,1 


2Pi/2,l,l 


25'i/2,l,l 


2^3/2,2,2 




V2(g+2) 
4 


^(9+2) „ 
4^3 







2P3/2,2,1 


4 ^+ 











2^3/2,1,1 


V2(g+2) 











2Pi/2,l,l 








9-4 „ 





2S'i/2,l,l 














-feo 



(Table 5a) 





2P3/2,2,0 


2P3/2,1,0 


2Pi/2,l,0 


2Si/2,l,0 


2Pi/2,0,0 


25i/2,0,0 


2^3/2,2,1 


V3(9+2) 
4 ^ 


4^3 













2^3/2,1,1 


_9±2 
12 ^ 


5(9 + 2) 
12 


6 





^(9-1)^ 
3 





2Pi/2,l,l 


V^(9-l)„ 
6 


v^(9-l) _ 
6 


9-4„ 










2Si/2,l,l 
















2^- 





2P3/2,2,1 


2P3/2,2,0 


V^{9+2) „ 
4 ^+ 


2P3/2,1,0 


4V^'^ + 


2Pi/2,l,0 




2Si/2,l,0 





2Pi/2,0,0 





25i/2,0,0 






(Table 5b) 

2P3/2,1,1 



2Pi/2,l,l 251/2, 1,1 



9+2 „ 
— — P ^ 
12 *= + 

5(9+2) 
12 '^+ 

^(9-1) 

6 "= + 



V^(9-l) „^ 
3 + 





V^(9-l) ■ 



6 

^(9-1) 
6 



9-4 „ 









^e+ 



(Table 5c) 



2P3/2,2,0 



2P3/2,1,0 



2Pi/2,l,0 



251/2,1,0 



2Pi/2,0,0 2Si/2,0,0 



2P3/2,2,0 

2P3/2,1,0 

2Pi/2,l,0 
2Si/2,l,0 
2Pi/2,0,0 
25'i/2,0,0 



V^(9-l) 



eo 



eo 



9+2 
6 





V2(g-1) 
3 





eo 



v^(9-i) 

3 






9-4 „ 



eo 



jeo 





v^(9-l) 
3 

9-4 



eo 



eo 







-feo 





(Table 5d) 
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2^3/2,2,-1 


2^3/2,1,-1 


2Pi/2,l,-l 


25i/2,l,-l 


2^3/2,2,0 

2^3/2,1,0 

2Pi/2,l,0 


V3(9+2) 


B+lp 

12 

5(9+2) 
12 

v^(9-l)„ 
6 


x/2(9-l) . 
6 

v^(9-l) _ 
6 







2Si/2,l,0 
2Pi/2,0,0 








v^(9-l)„ 
3 




9^e 





25i/2,0,0 














(Table 5e) 





2^3/2,2,0 2P3/2,1,0 2Pi/2,l,0 2Si/2,l,0 2Pi/2,0,0 2Si/2,0,0 


2^3/2,2,-1 
2^3/2,1,-1 

2Pi/2,l,-l 

251/2,1,-1 


V^(9+2) f^^Ci 

«+/e+ '^'i2'^e+ ^(^^'6+ ^(f^)e+ 
^<r''e+ ^'r''e+ «-^e+ 
fe+ fe+ 



(Table 5f) 





2P3/2,2,-l 2P3/2,1,-1 2Pi/2,l,-l 2Si/2,l,-l 


2P3/2,2,-2 


2^3/2,2,-1 
2^3/2,1,-1 

2A/2,1,-1 
2Si/2,l,-l 


-f^'^o -f^eo 

-^-0 -§75^" -"t^^" 
feo 


V^(9+2) „ 
V^(9+2) „ 




2P3/2,2,-2 


v^(9+2) v^(9 + 2) S-lr^ 





(Table 5g) 
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